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In this research paper we will cover the study of the modelling and the 
control methods of the variable speed wind turbine based on doubly-fed 
induction generator (DFIG). It represents the most stressed structure given 
its distinctive characteristics. To control the electrical powers generated by 
this system independently, the vector control with the stator flux orientation 
is founded according to two techniques: i) the control of the powers by the 
backstepping technique and ii) the robust control based on the active 
disturbances rejection control. After the synthesis of the controllers of those 
two methods, their performances will be tested and compared to evaluate 
their effectiveness. We are mainly interested in the robustness test of the two 
control strategies with respect to the internal parameters’ fluctuation of the 
generator. The computing of the different parameters regulators of these two 
strategies is carried out using a genetic algorithm. This computing method 
makes it possible to arrive at an optimal solution of the DFIG power control. 


Wind energy The different parts are simulated using MATLAB/Simulink environment. 
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NOMENCLATURE 

R : Turbine radius 

P : Air density 

v : Wind speed 

Qi : Turbine speed 
Cp : Power coefficient 
À : Velocity ratio 

B : Pich angle 

J : Total inertia 
Jrurbine : Turbine inertia 
Jo : Generator inertia 


f : Coefficient of viscous friction 
G : Gain multiplier 
saq : dq axis Stator fluxes 


raq : dq axis Rotor fluxes 

Vsaq  : Stator voltages according to dq axis 
Vaq  : Rotor voltages according to dq axis 
İsaq  : Stator currents according to dq axis 
iraq __: Rotor currents according to dq axis 
Ws : Rotor angular frequency 

wr : Stator angular frequency 


Journal homepage: http://ijpeds.iaescore.com 
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1. INTRODUCTION 

Faced with the ever-increasing demand for energy, industrialized countries have made massive use 
of nuclear power plants. This type of energy has the advantage of not causing atmospheric pollution, but the 
risk of nuclear accident [1], [2]. Nowadays countries are increasingly turning to the use of renewable energy 
sources as a solution to that problem. Wind power in term of example, which represent an attractive 
alternative that may in the future replace traditional power plants [1], [3], [4]. 

The variable speed wind turbine based on a doubly-fed induction generator (DFIG), driven by a 
turbine and controlled by the rotor is the most used system in wind energy conversion considering the 
advantages it offers [3], [5]-[7]. In this structure, the stator of the generator is connected directly to the grid, on 
the other hand the rotor is linked to the grid via a transformer and link filter with two bidirectional 
converters [3], [4], [8], [9]. Figure 1 shows the functional scheme of this structure based on the DFIG and two 
power converters. 

For this system command, various approaches are available in the literature. A study by Nadour and 
Essadki [10] shows that the backstepping approach provides good results compared to the proportional integral 
(PI) control especially with regard to the hounding the references values and robustness against external and 
internal disturbances. However, this approach does not present more significant results faced with the change 
of internal DFIG parameters with higher variation rates (100% for example). Especially since it does not 
present an analytical method allowing the computing of the regulator parameters. 


Wind turbine Utility Grid 
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Figure 1. Wind turbine structure based on DFIG [1] 


Bossoufi and Karim [11] have presented a control technique based on a control approach of pole 
placement incorporated with a backstepping control system. The experimental results obtained based on a 
prototyping platform realized with DFIG-generator, field programmable gate array (FPGA) and wind turbine 
[11] show the efficiency of the technique for the system stability. The problem of this solution is that it does not 
present any detail about the robustness in the case of internal parameters' variation. Indeed, these parameters 
may undergo changes caused by many physical phenomena (temperature, inductors saturation) [12]. 

Rouabhi and Rachid [13] presented a combination of vector and backstepping control. Simulation 
results show the decoupled control of both electrical powers and a good performance likened to the results 
acquired in the case of the compensator use. However, we do not find in this approach a method for the 
synthesis of the controllers’ parameters, on which the performances of this approach are based. 

To ensure system stability and rapid dynamic response, the control parameters must be well chosen. 
However, in all the previous works and in the majority of the documents which treat the same subject, one 
does not find a study which presents a procedure for the determination of these control strategy parameters. 
The calculation is often done by a random adjustment during the simulation. 

In recent years, several works have dealt with the DFIG command used in wind energy conversion 
system (WECS) by active disturbance rejection control (ADRC) [14]-[16]. With this control method, the 
system is able to compensate external disturbances and those due to parameter changes. The work of this 
research paper is focused on the study and comparison of this methods robustness with that based on 
Backstepping. We also propose a method based on genetic algorithms, for the calculation of the different 
controllers’? parameters, which allows to arrive at an optimal solution of DFIG power control. We can show 
that replacing backstepping controllers with ADRC controllers optimized by genetic algorithm (GA) allows to 
improve significantly the system performances. 
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This paper is structured as shown in: the first part is devoted to the studied system modeling with its 
different components. In the next part, we treat the DFIG control principle by the backstepping technique. The 
third part is entirely devoted to the presentation of the DFIG robust control founded on the ADRC. The next 
part is devoted to the presentation of the procedure used for the optimization of the two control approaches 
using the genetic algorithm. The last section of this paper is devoted to the presentation and analysis of the 
simulation results, and the robustness test using the MATLAB/Simulink environment. 


2. RESEARCH METHOD 
2.1. Modelling of wind turbine using DFIG 
2.1.1. Wind turbine model 
The wind turbine converts the kinetic wind energy and drives the variable speed generator. The 
maximum power which can be obtained according to Betz's law is [17]-[19]: 


1 
Pextr =5p mR? v°C,(B,A) (1) 


where C, is the power coefficient, which is defined by the ratio between the total extracted power and the 
incident one [17]. This coefficient depends on A, the velocity ratio, and f the pich angle, it is formulated by the 
following function [16], [20]: 


0.035 


1 1 0.035 -21(-—+—+ ) 
C(b) = 5 (1166 Page Oe 5) e  4+0.088 1+6% + 0.00684 (2) 


the aerodynamic torque Taer is given by: 


Pextr 1 
= A = m? R? v3C,(B,A) (3) 
dQ.mec 
JT- + fAmec = Ty = Tem (4) 
J = Turbine +Jg (5) 


from (1), (3), and (5), we can deduce the functional diagram of the turbine model which is illustrated in this 
Figure 2. 


Qmec 


Figure 2. Functional scheme of turbine model 


2.1.2. DFIG modelling 
The equations which describe the electric part of the generator in the two-phase reference are given 
by the following relations [12], [21]: 


Psa = Lslsa + Lmlra 
Psq = Lslsq + Linlrq 


6 
Pra = Lyla + Lmlsa (6) 
Prq = Lylrg + Lmlsq 
dos 
Vsa = Rslsa + “Esd — Ws Psq 
dos 
Voq = Rslsq + a + WsPsa 7) 


_ dfrd = 
Vra = Rylra T dt Wr Prq 


do 


Vaq = Rrlrg +7 + @rPra 
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the equation describing the mechanical part is as shown in [21]: 

Tem = PÈ (Ira Psq ~ Inq Psa (8) 
where p=2 is the number of pole pairs. We have considered the stator flux control oriented along the d axis of 


the two-phase reference. It is also considered that the stator resistance Rs is negligible during the model 
computing [12], [17], [21], [22]. So, we deduce: 


Psa = Ps = cte 
f Psq = 0 (9) 
Vsa 
a = Ws Ps 
dor 
Vig = Rrra + i = Wr Prq (10) 
dor 
Veq = R,Irq + ona + Wr Pra 
Psa = Ps = Lslsa + Lmlra 
0 = Lsls + Lmlrq an 
Pra = Lrlra + Lmlsa 
Prq = Lrlrq + Lmlsq 
S Lm 
Isa = - T q, ra 
_ _ Lm (12) 
I5q = Tla 
Tem = -+p I (13) 


Finally, we can deduce the block diagram of the DFIG model (Figure 3). 
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Figure 3. Block diagram of the DFIG model 


2.2. Powers control using the backstepping technique 
2.2.1. Backstepping principle 

The principle of backstepping relies on its ability to dissociate a complex nonlinear problem control 
for a system into simpler ones. Basically, this technique synthesis is separated into different design steps. At 
each one, we're virtually dealing with a simpler system with single-input and single-output problem, and each 
sub system supply a baseline for the next synthesis step. The global stability and overall performance are 
obtained by a Lyapunov function for the global system [4], [6], [11], [13]. 


2.2.2. DFIG control based on the backstepping method 

The backstepping method that we are going to apply for controlling the double-fed asynchronous 
machine is based on the vector control. In this approach, the command will be established in the reference (d 
and q) according the condition in (11). The control of the rotor side converter (RSC) allows the regulation of 
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the exchanged powers with the grid. the synthesis of the control blocks with the backstepping method is 
divided into two subsystems: the first one is to generate the current reference I,a and Iq and the second one is 
to generate the voltage reference V,a and Vu [4], [10], [23]. the relations linking the stator powers and the 
rotor currents are expressed by (14) and (15): 


pref = Vslm pref 


S Ls rq 
me ee (14) 
S = Ws.Ls Ls ` rd 
E R. 
sm (15) 
rer _ Ls ref 
rd Vg. Lm “S 
the rotor currents derivatives can be mathematically determined and given by the following: 
dI 1 
ae = (Vya — Rr- lra + 9. Ws. Lr- 0. Lados (16) 
dl LmVs\ 1 
8 = (Vg — Rr- Irq — G- Ws. by. O. Lra — 9.Ws- aa ae (17) 


L? ; $ ates F es 
where o = (1 — AE) represent the dispersion coefficient. The computing of the reference quantities of the 
SST 


currents and voltages control is based on the Lyapunov functions according to two steps [4], [11]-[13]. 
- Step 1: the errors E; and E2 between the currents Irq, Ira and their reference Trgrer and Trdret is given by (18). 


E = eet lr 


(18) 
E, = ig — hra 
The derivative of this error is given by (19). 
SN ref ` 
A= e -ha (19) 
; re ; 
Ez = La — ird 
We select the Lyapunov function so such that: 
v = > (E? +E?) (20) 
The errors are chosen so as to have a zero derivative of the Lyapunov function: 
E, = —K,E, where: K, > 0 (21) 
E, = —K,E, where: K, > 0 (22) 
v = —K,E? — KE? (23) 
. Le A 1 Lim Vs 
È = (=e?) ae eq — Rr drg — 9-We- by Obra — 9-8) (24) 
: Ls : 1 
È, = Ga z7) - Tg Vea Reha + g-Ws- by. O-Irq) (25) 
- Step 2: 
Ls 2 1 1 LmVs 
—K,E, = Gee) = Veg — og Reng — Wee by Oleg — 9.8) (26) 
ls 2 1 1 
—K,E> = Gar or) = Vea — Ao (Re bra + g Wo. br 0. Irq) (27) 


From the previous equations, the voltages reference can be mathematically determined: 
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La `i LmVs 
Vig = Ly Gata + K,E,) + Rr- lq + g-Ws- Lr. Olya + 9. E (28) 
Ls ; 
Vya = Ly.o (- A OF! +K2E3) + Rp- lya + J-Ws- Lr- 0- lrq (29) 


Kı and Ko are the backstepping controller parameters, they are strictly positive, and ensure the system 
stability as well as a fast dynamic response [14]. These performances are ensured with a zero derivative of 
Lyaponov function. The functional scheme of the backstepping command is given in the Figure 4. 
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Figure 4. Block diagram of the control structure by backstepping 


2.3. Active disturbance rejection control strategy 
2.3.1. Active disturbances rejection control strategy 

The ADRC is a robust control strategy which is based on the model system extension using a state 
observer extended state observer (E.S.O) to have an estimation and to cancels, in real time, all disturbances 
affecting the system whether internal or external [15], [16], [18], [24]-[26]. To explain the principle of the 
ADRC, we treat the case of a first order system: 

of = —<y(t) + b.u(t) where: b = by + Ab and by = = (30) 
where b is a parameter to be estimated, bo represents the known part of b and Ab represents the modeling 
error due to the system parameters changes. u and y are input and output variables, K is the static gain and T 
is a system constant. 


oo = -5y + Ab.u(t) + bo. u(t) 


dy(t) _ 
a f(y, d,t) + bo.u(t) 


where: f(y, d,t) = —<y(t) + =d(t) + Ab. u(t) = y(t) (31) 
f(y,d,t) represents the effect resulting from internal parameters variations and the external disturbances 
impact. The basic idea is the estimation and compensation of the disturbance. The (31) can be written in an 


augmented state space form as [27]: 


Xy = X, +bo.u 


%=h  whereh=f (32) 
y = xX 
x = Ax + bBu + Eh 


where: A = G o) B = (*), C = (10) and E = (D. The observer will be able to estimate the derivatives of y 
and f and it is developed as shown in: 


S ae (34) 
ğ = Cx 
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where L represent the observer’s gain vector. 


In Figure 5(a), we represent the structure of the ADRC controller previously studied [14], [15], and the 
structure of the extended state observer is given by the Figure 5(b) [28]. 


|a bou 
Uo 1 u | 
r —>Q— Kp <Q i Process : 
p 5 y e =R Ry 
| —>| 2 i 
s 


X2 


ule 


Figure 5. Structure of (a) the ADRC controller and (b) ESO 


2.3.2. Control of the rotor side converter 

By referring to the assumptions relating to the stator flux orientation along the d-axis, and neglecting 
the stator resistance, the equations of the mathematical model of DFIG can be simplified [16]. So, the rotor 
currents are expressed by the following relations: 


dird Rr . 2 1 

= —— ig + wr lrg +V, 35 
ge gh me Or ira tog, Vra (35) 
dirq Ry . . 1 
JF tog — Ore lra — Or aT —V, 
dt oL, 14 rra T oLyLs S oLy rq (36) 


we can also put these relations in the (37) and (38). 


Ry . Fi 1 
= ——irg + Gy. lrg + a bo)Vra 


; d 
d ; L 
E = fa(i,a d, t) + bou (t) where: an i i i (37) 
u = Vra , Do = 
oLr 
PE Rr; = = Lm ae _— 
dirg _ T oLr ‘rq Or: lra r oLrLs Ps i Ce borg 


rn pees d, t) + b ou (t) where: 1 (38) 
u= Veg , bro = OL, 
fy and f, represent the disturbances impacting the rotor currents Iq and Iq. Via and Vq are respectively the 
command variables of the current’s loops Ira and Irq. bro represent the known part of the process parameters. 
The schematic diagram of the rotor side control is given by Figure 6. The parameters Ky, Bı and B2 are 
determined to have a perfect follow-up of Iar and Iq references values. These parameters are obtained in our 


case using GA. 


2.3.2. Control of the GSC 
The currents igr and iaf are given by (39). 


digf _ 1 Rf. P 1 
dt = L Vea = ie laf = Wslqf = rm (39) 
digf _ 1 Rf. r 1 
ae gy ee g a Wala ip Var (40) 


The previous relations can be put in the (41) and (42). 


= Rf. ; 1 
diaf far =g Sapar Wst Cam Na 
= fap (iap 4, t) + bou(t) where: 


(41) 


dt 1 
u = Var, bro =- 
af , Pfo P 
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_ 1 Rf. n 1 
far = Vsa i a slag * G = bo War 


di f Ly 
a = fas apes t) + bou(t) where: 


i (42) 
Vea ARa 


The computing of the parameters used in the regulation algorithm, is done to have a perfect tracking of the reference 
values by the currents iqf and igs. This control algorithm by ADRC can be summarized according to Figure 7. 


larref Havre 
I 


dr 


I qrRef 


6s ——— PLL |<— ¥s 


Figure 6. The GSC control diagram using ADRC Figure 7. The GSC control diagram using ADRC 


2.4. Genetic algorithm 

The genetic algorithm is an optimization method for finding optimal solutions to research problems. 
It is based on the natural selection methods. This approach was proposed by John Holland in 1975, and then 
adopted as a very efficient method to find the solutions to an optimization problem. It makes it possible to 
avoid the local minima constituting a major problem in the case of nonlinear systems [29], [30]. 

The basic principle is based on the repeated modification of the population so as to evolve towards 
an optimal solution, and this starting from an initial population consisting of a certain number of elements 
appointed chromosomes. The relevance of each solution is measured by a fitness function. The algorithm is 
based primarily on three stages: selection, crossing, and mutation. These three operations must be applied to 
create new individuals to be evaluated with respect to their parents. The algorithm stops when we obtain 
individuals with optimal performances [29], [31], [32]. 

In the case of DFIG powers’ control, the objective function can be computed and evaluated for each 
individual belonging to the actual population. Therefore, the parameters of the ADRC or Backstepping 
controller can be decoded for each individual, then the DFIG model is simulated to obtain the objective 
function value which consists in following the reference power required by the turbine. The following 
population is generated based on the genetic algorithm operators (selection, crossing and mutation) [30], [33]. 
These two steps are repeated from generation to generation until a final population is obtained. In Figure 8, 
we present the optimization diagram to be carried out to reach the optimal parameters values of the 
backstepping and ADRC regulators. The types of the GA operations used and the parameters’ initialization are 
given in Tables 1 and 2. In this work, the power of the machine used is 1.5 MW. The algorithm stops when the 
best fitness value is less than or equal to the limit fitness value which has been set here at 5.104 MW. 
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Figure 8. Optimization block diagram based on genetic algorithm 
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Table 1. Types of the GA operations used 


Property Type 

Selection GA default selection function 
Mutation Uniform 
Crossover Arithmetic 


Table 2. Genetic algorithm parameters 


Property valie 
Backstepping _ADRC 
Variables to optimize 2 3 
Population size 50 50 
Maximum size of generations 500 500 
Mutation fraction 0.01 0.01 
Crossover fraction 0.08 0.08 
Tolerance 5.10? 5.10? 


2.4.1. Parameter’s optimisation using GA for backstepping technique 

As mentioned before, the control system must be able to guarantee the two main objectives: closed 
loop stability of the system and zero regulation of the tracking error. The problem which arises at this level is 
the determination of the backstepping adjustment coefficients Kı and K2. To avoid determining the gains of 
the regulators with non-analytical methods (random), we propose in this work an optimization by the genetic 
algorithm. The principle consists in considering an adaptation law established by the Lyapunov function as a 
minimization problem. This function is given as presented before by (43) [10], [11], [13]. 


py = —K,E? -RE (43) 
Where E = i — Iq and E = S — l,a. In fact, the control of the I and Ir currents allows the 


generation of an active power P in accordance with the reference power required by the mechanical part 
(turbine). A GA based backstepping control optimization can be represented as shown in Figure 9. 


Prof 


' ' 
P mes le 
Figure 9. GA-based backstepping control optimization 


2.4.2. Parameter’s optimisation using GA for ADRC 

The ADRC allows an estimation, using a state observer, of modeling uncertainties and external 
disturbances, and a compensation for them in real time [20]. The objective is to have a faster observer 
dynamic than the closed loop system one, as well as a reduced noise sensitivity [15], [16], [18], [28]. The 
genetic algorithm is used to determine the Ky, Bı and B2 parameters. These parameters are adjusted so as to 
obtain the optimum performance in terms of dynamics, robustness and disturbance’s rejection. The objective 
function is therefore chosen to minimize the difference between the measured power and the reference one. A 
GA based ADRC optimization can be represented in Figure 10 and the parameters values used for genetic 
algorithm in this work are presented in Table 2. 


ref _ _ Vs.L ref 

pret = — Yelm re (4) 
ref _ vè _Yslm yref 

D =e p (45) 
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Figure 10. GA-based ADRC optimization 


The best score value and mean score versus generation for the two control methods are represented 
in Figure 11. The parameters’ evolution trought the generation is presented in Figure 12. The backstepping 
gains K; and Ko, and the ARDC controller parameters Kp, B; and Bz are listed in Table 3. 
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Figure 11. Backstepping and ADRC fitness value 
Table 3. Backstepping and ADRC parameters’ values 
arameter mbo. alue obtaine ower range er range 
Paramet Symbol Value obtained by GA L U 
Backstepping gain 1 Kı 38824 1.0 50000 
ackstepping gain 2 g 
Backst 2 K: 29986 1.0 50000 
Currents controller gain (ADRC) K, 17175 1.0 20000 
xtended state observer gain 1 . : 
Extended state ob: s (ADRC 4357.5 1.0 5000 
B2 481940 1.0 500000 
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Figure 12. Parameters’ evolution through the generation 
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3. RESULTS AND DISCUSSION 
3.1. Tracking and control test 

The parameters of the DFIG are listed in Table 4. In Figure 13 we represent the wind profile used in 
this work which is modelled by (46): 


v(t) = 8 + 2sin (2.5¢ ? z) + 2sin (4t = z) + 1.5sin (5.4t = =) + 0.5sin (2.5¢ 2 Z) (46) 

Figures 14 and 15 show the simulation results of the produced powers and the rotor currents Ira and 
I,q with their references, for the Backstepping technique and ADRC. The reference reactive power is chosen 
zero to guarantee a unit power factor. As can be seen in the previous figures, the produced powers controlled 
by backstepping and ADRC, optimized by genetic algorithm, follow perfectly their references. We can notice 
more precision and stability with a response without fluctuations or exhibiting overshoot on that based on the 
backstepping technique. From these figures, we can see that the stator currents vary according to the wind 
profile provided for the simulation and the frequency of the rotor voltages follows the DFIG speed variation. 
In Figure 16, we present the stator currents and rotor voltages of the DFIG. 


Table 4. List of DFIG parameters 


Parameter Symbol Value 
Rated power Ps 1.5 MW 
Stator resistance Rs 0.0089 Q 
Rotor resistance Rr 0.0137 Q 
Stator inductance Ès 0.0137 H 
Rotor inductance E 0.01367 H 
Mutual cyclic inductance Ei 0.0135 H 
Number of pole pairs p 2 
Optimal tip speed ratio of turbine opt 8.1 
Maximal power coefficient Comax 0.48 
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Figure 13. The wind speed profile 
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Figure 14. Active and reactive power 
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Figure 16. Stator current and rotor voltage 


3.2. Test of the robustness 

The objective of the robustness test is to verify the control methods performance against the internal 
parameters variation of the DFIG. In fact, the regulator parameters computing is based on a DFIG model 
whose parameters are considered to be fixed. We are more particularly interested in highlighting the 
performance limits of the backstepping method in this case. However, these parameters can change over time 
or because of external phenomena. In Figures 17 and 18 we represent the produced powers evolution and the 
rotor currents I,q and Ia after a variation R, and L, with a rate of 100%. 

From the Figures 17 and 18, the variations of the R, and L, values did not have much impact on the 
performance of the system controlled by the active disturbances rejection control unlike the backstepping 
technique. Indeed, the ADRC allows the estimation and compensation of disturbances. The tracking of the 
reference values is always assured as well as the stability which is not affected. 
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Figure 17. Produced powers after a 100% variation of Rr and Lr 
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Figure 18. Currents Irq, Ira and its reference after a 100% variation of Rr and Lr 


4. CONCLUSION 

In this research article, we have dealt the modeling and control of the wind turbine driven by DFIG. 
We are interested to highlight and evaluate the dynamic performances of ADRC and backstepping methods 
and their sensitivity to the DFIG’s internal parameters variation. The two control methods studied consist in 
controlling the rotor currents in order to adapt the DFIG rotation speed to the wind speed, by the control of 
the electromagnetic torque, and also in controlling the active and reactive powers injected into the grid. The 
synthesis of the regulation parameters is based on the genetic algorithm which allows to have the optimal 
values to guarantee the best performances. According to the simulation results, and for normal operation 
assuming that the DFIG's internal parameters are fixed, we notice that the backstepping control, optimized by 
GA, allows good performances, in particular a fast response and without overshoot, unlike the ADRC which 
has some undulations at the start. However, internal parameters are always subject to variations. In that case, 
we have found that the ADRC allows a good tracking of the desired values as well as a robustness against the 
DFIG’s internal parameters variation. In terms of perspectives, the rate of the internal parameter’s variation 
treated in this work is equal to 100%. However, we do not find any work that studies the order of magnitude 
of this variation rate, in order to have more precise values. We are therefore considering the development of 
methods or algorithms making it possible to estimate this variation rate and take it into account in the control 
laws used. 
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